Wormhole geometries in f{R) modified theories of gravity 
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In this work, we construct traversable wormhole geometries in the context of f{R) modified the- 
ories of gravity. We impose that the matter threading the wormhole satisfies the energy conditions, 
so that it is the effective stress-energy tensor containing higher order curvature derivatives that 
is responsible for the null energy condition violation. Thus, the higher order curvature terms, in- 
terpreted as a gravitational fluid, sustain these non-standard wormhole geometries, fundamentally 
different from their counterparts in general relativity. In particular, by considering specific shape 
functions and several equations of state, exact solutions for f{R) are found. 
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I. INTRODUCTION 

Various independent high-precision observational data 
have confirmed with startling evidence that the Universe 
is undergoing a phase of accelerated expansion [l| . Sev- 
eral candidates have been proposed in the literature to 
explain this phenomenon, ranging from dark energy mod- 
els to modified theories of gravity. In the latter context, 
one may assume that at large scales Einstein's theory of 
general relativity breaks down, and a more general action 
describes the gravitational field. The Einstein field equa- 
tion of general relativity was first derived from an action 
principle by Hilbert, by adopting a linear function of the 
scalar curvature, i?, in the gravitational Lagrangian den- 
sity. However, there are no a priori reasons to restrict 
the gravitational Lagrangian to this form, and indeed 
several generalizations have been proposed. In particu- 
lar, a more general modification of the Einstein-Hilbert 
gravitational Lagrangian density involving an arbitrary 
function of the scalar invariant, f{R), was considered in 
0, and further developed in 

In this context, a renaissance of f{R) modified theories 
of gravity has been verified in an attempt to explain the 
late-time accelerated expansion of the Universe (see Refs. 
[J for a review). Earlier interest in f{R) theories was 
motivated by inflationary scenarios as for instance, in 
the Starobinsky model, where f{R) = R — A + aR? was 
considered In fact, it was shown that the late-time 
cosmic acceleration can be indeed explained within the 
context of f{R) gravity [6]. Furthermore, the conditions 
of viable cosmological models have been derived Q , and 
an explicit coupling of an arbitrary function of R with 
the matter Lagrangian density has also been explored 



[a] . Relative to the Solar System regime, severe weak field 
constraints seem to rule out most of the models proposed 
so far [13] , although viable models do exist [ll[ . In the 
context of dark matter, the possibility that the galactic 
dynamics of massive test particles may be understood 
without the need for dark matter was also considered in 
the framework of f{R) gravity models [T^ . 

The metric formalism is usually considered in the lit- 
erature, which consists in varying the action with respect 
to g^'^ . However, other alternative approaches have been 
considered in the literature, namely, the Palatini formal- 
ism [l^ where the metric and the connections are 
treated as separate variables; and the metric-affine for- 
malism, where the matter part of the action now depends 
and is varied with respect to the connection • The ac- 
tion for /(i?) modified theories of gravity is given by 

-g f{R) + SM{9^'',^), (1) 

J 

where k, = SvrG; throughout this work we consider k = 1 
for notational simplicity. SM{g^'^ is the matter ac- 
tion, defined as Sm = J d'^x^J—g Lmig , 'il^) , where Lm 
is the matter Lagrangian density, in which matter is min- 
imally coupled to the metric g^i, and ?/; collectively de- 
notes the matter fields. 

Now, using the metric approach, by varying the action 
with respect to g^^ , provides the following field equation 

FR^.. ~ \fgt.. - V^V.F -f g^.UF = T™ , (2) 

where F = df/dR. Considering the contraction of Eq. 
(121), provides the following relationship 



FR-2f + iUF = T, 



(3) 
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which shows that the Ricci scalar is a fully dynamical 
degree of freedom, and T = T^f^ is the trace of the stress- 
energy tensor. 



2 



In this work, we extend the analysis of static and 
spherically symmetric spacetimes considered in the lit- 
erature (for instance, see [11]), and analyze traversable 
wormhole geometries in /(i?) modified theories of gravity. 
Wormholes are hypothetical tunnels in spacetimc, possi- 
bly through which observers may freely traverse. How- 
ever, it is important to emphasize that these solutions are 
primarily useful as "gedanken-experiments" and as a the- 
oretician's probe of the foundations of general relativity. 
In classical general relativity, wormholes are supported 
by exotic matter, which involves a stress-energy tensor 
that violates the null energy condition (NEC) [l6l [13] ■ 
Note that the NEC is given by T^,i,k>'k'' > 0, where fc'' is 
any null vector. Thus, it is an important and intriguing 
challenge in wormhole physics to find a realistic matter 
source that will support these exotic spacetimes. Several 
candidates have been proposed in the literature, amongst 
which we refer to solutions in higher dimensions, for in- 
stance in Einstein-Gauss-Bonnct theory [18, 19], worm- 
holes on the brane [23]; solutions in Brans-Dicke theory 
[2ll . [2^ . [23| : wormhole solutions in semi-classical gravity 
(see Ref . [2J| and references therein) ; exact wormhole so- 
lutions using a more systematic geometric approach were 
found [2^; geometries supported by equations of state 
responsible for the cosmic acceleration [2^, solutions in 
conformal Weyl gravity were found [27*] , and thin accre- 
tion disk observational signatures were also explored [2^ , 
etc (see Refs. [2^, [S^] for more details and [13] for a re- 
cent review). 

Thus, we explore the possibility that wormholes be 
supported by f{R) modified theories of gravity. It is 
an effective stress energy, which may be interpreted as a 
gravitational fluid, that is responsible for the null energy 
condition violation, thus supporting these non-standard 
wormhole geometries, fundamentally different from their 
counterparts in general relativity. We also impose that 
the matter threading the wormhole satisfies the energy 
conditions. 

This paper is organized in the following manner: In 
Sec. ini the spacetime metric, the effective field equations 
and the energy condition violations in the context of f{R) 
modified theories of gravity are analyzed in detail. In 
Sec. mil specific solutions are explored, and we conclude 
in Sec. HVl 



II. WORMHOLE GEOMETRIES IN f{R) 
GRAVITY 



where $(r) and b{r) are arbitrary functions of the radial 
coordinate, r, denoted as the redshift function, and the 
shape function, respectively [l^. The radial coordinate 
r is non-monotonic in that it decreases from infinity to a 
minimum value tq, representing the location of the throat 
of the wormhole, where b{ro) — tq, and then it increases 
from ro back to infinity. 

A fundamental property of a wormhole is that a flaring 
out condition of the throat, given by {b — b'r)/b'^ > 0, is 
imposed [IB], and at the throat 6(ro) — r = tq, the con- 
dition b'{ro) < 1 is imposed to have wormhole solutions. 
It is precisely these restrictions that impose the NEC vi- 
olation in classical general relativity. Another condition 
that needs to be satisfied is 1 — b{r) jr > 0. For the worm- 
hole to be traversable, one must demand that there are 
no horizons present, which are identified as the surfaces 
with e'^* 0, so that ^{t) must be finite everywhere. 
In the analysis outlined below, we consider that the red- 
shift function is constant, = 0, which simplifies the 
calculations considerably, and provide interesting exact 
wormhole solutions (If ^ 0, the field equations be- 
come forth order differential equations, and become quite 
intractable). 

The trace equation ^ can be used to simplify the field 
equations and then can be kept as a constraint equa- 
tion. Thus, substituting the trace equation into Eq. ([2]), 
and re-organizing the terms we end up with the following 
gravitational field equation 
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(5) 



where the effective stress-energy tensor is given by — 



rlfj + fjZ'\ The term fjZ'^ is given by 



(6) 



(c) 

and the curvature stress-energy tensor, , is defined as 



1 



1 



V^V.F - -g^, {RF + aF + T) 



(7) 



It is also interesting to consider the conservation law 
for the above curvature stress-energy tensor. Taking into 
account the Bianchi identities, V^G^jy = 0, and the dif- 
feomorphism invariance of the matter part of the action, 
which yields VT^™' = 0, we verify that the effective 
Einstein field equation provides the following conserva- 
tion law 



A. Spacetime metric and gravitational field 
equations 



(8) 



Consider the wormhole geometry given by the follow- 
ing static and spherically symmetric metric 



1 -b{r)/r 



+r^(dr+sin^6'd0^), (4) 



Relative to the matter content of the wormhole, we im- 
pose that the stress-energy tensor that threads the worm- 
hole satisfies the energy conditions, and is given by the 
following anisotropic distribution of matter 



T^,i, = (p + Pt)Uf, Uu + Pt gtj.i' + {Pr 



Pt)Xt^Xl': (9) 
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where [/^ is the four- velocity, is the unit space hke 
vector in the radial direction, i.e., = -y^l — b{r)/r S^r- 
p{r) is the energy density, Pr{r) is the radial pressure 
measured in the direction of x^ i s-nd pt(r) is the trans- 
verse pressure measured in the orthogonal direction to 
X''. Taking into account the above considerations, the 
stress-energy tensor is given by the following profile: 
T^'^ = d\a.g[- p{r),pr{r),pt{r),pt{r)]. 

Thus, the effective field equation ^ provides the fol- 
lowing relationships 



(10) 
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= ^^ 
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F" - F' 
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b'r-b 
2r2(l - b/r)_ 



- H 













r 



(11) 



(12) 



where the prime denotes a derivative with respect to the 
radial coordinate, r. The term H — H{r) is defined as 



H{r) ^- (FR + UF + T), 



for notational simplicity, 
given by 

R 



The curvature scalar, 
26^ 



and DF is provided by the following expression 



UF = 



r 



F" - 



2r2(l - b/r) 



, 2F' 
F' + 

r 



(13) 
i?, is 

(14) 
(15) 



Note that the gravitational field equations PTH) - P^ . 
can be reorganized to yield the following relationships: 



Fb' 



Pt = 



bF_ 

F' 

r 



^{b'r-b)-F" 



b 

1 - - 
r 



b 

1 - - 

r 



F 

2^(^ 



b'r).. 



(16) 
(17) 

(18) 



which are the generic expressions of the matter thread- 
ing the wormhole, as a function of the shape function 
and the specific form of F{r). Thus, by specifying the 
above functions, one deduces the matter content of the 
wormhole. 

One may now adopt several strategies to solve the field 
equations. For instance, if b{r) is specified, and using a 
specific equation of state Pr = Pr{p) or pt = Pt{p) one can 
obtain F{r) from the gravitational field equations and 
the curvature scalar in a parametric form, R{r), from its 
definition via the metric. Then, once T = ^ is known 
as a function of r, one may in principle obtain f{R) as a 
function of R from Eq. ([3]). 



B. Energy condition violations 

A fundamental point in wormhole physics is the en- 
ergy condition violations, as mentioned above. However, 
a subtle issue needs to be pointed out in modified theories 
of gravity, where the gravitational field equations differ 
from the classical relativistic Einstein equations. More 
specifically, we emphasize that the energy conditions 
arise when one refers back to the Raychaudhuri equa- 
tion for the expansion where a term R^^k^k'^ appears, 
with fc^ any null vector. The positivity of this quantity 
ensures that geodesic congruences focus within a finite 
value of the parameter labelling points on the geodesies. 
However, in general relativity, through the Einstein field 
equation one can write the above condition in terms of 
the stress-energy tensor given by T^^k^k'^ > 0. In any 
other theory of gravity, one would require to know how 
one can replace R^i, using the corresponding field equa- 
tions and hence using matter stresses. In particular, in a 
theory where we still have an Einstein-Hilbert term, the 
task of evaluating R^^k^k'^ is trivial. However, in f{R) 
modified theories of gravity under consideration, things 
are not so straightforward. 

Now the positivity condition, Rf^^k^k'^ > 0, in the 
Raychaudhuri equation provides the following form for 
the null energy condition T^^k^k" > 0, through the mod- 
ified gravitational field equation ^ , and it this relation- 
ship that will be used throughout this work. For this 
case, in principle, one may impose that the matter stress- 
energy tensor satisfies the energy conditions and the re- 
spective violations arise from the higher derivative cur- 

(c) 

vature terms Tf^J . Another approach to the energy con- 
ditions considers in taking the condition T^^k^k'^ > at 
face value. Note that this is useful as using local Lorentz 
transformations it is possible to show that the above con- 
dition implies that the energy density is positive in all lo- 
cal frames of reference. However, if the theory of gravity 
is chosen to be non-Einsteinian, then the assumption of 
the above condition does not necessarily imply focusing 
of geodesies. The focusing criterion is different and will 
follow from the nature of R^^k^^k'^ . 

Thus, considering a radial null vector, the violation of 
the NEC, i.e., T'^^ kt'k" < takes the following form 



P'^'+Pf 



F 



1 



F" - F'- 



b'r 



2r2(l - b/r) 

(19) 

where p'^^ + < 0. Using the gravitational field equa- 
tions, inequality (jl9p takes the familiar form 



-pf = 



b'r-b 



(20) 



which is negative by taking into account the flaring out 
condition, i.e., {b'r — b)/b^ < 0, considered above. 
At the throat, one has the following relationship 



P 



off 



cff| _ P + Pr 
Pr \ro — TP 



1 - b'iro) F' 



2ro 



F 



<0. (21) 
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It is now possible to find the following generic rela- 
tionships for F and F' at the throat: Fq < —2ro{p + 

Pr)|ro/(l-&') if > 0; andF^ > -2ro(p+p,.)|ro/(l " 
if F < 0. 

Consider that the matter threading the wormhole 
obeys the energy conditions. To this effect, imposing 
the weak energy condition (WEC), given by p > and 
P + Pr > 0, then Eqs. P^ -pT )) yield the following in- 
equalities: 



{2F + rF'){b'r -b) 
2^ 



F" 1 



Fb' 
b 



> 0. 



>0, 



(22) 
(23) 



respectively. 

Thus, if one imposes that the matter threading the 
wormhole satisfies the energy conditions, we emphasize 
that it is the higher derivative curvature terms that sus- 
tain the wormhole geometries. Thus, in finding wormhole 
solutions it is fundamental that the functions f{R) obey 
inequalities (dH) and 



III. SPECIFIC SOLUTIONS 

In this section, we are mainly interested in adopting 
the strategy of specifying the shape function b{r), which 
yields the curvature scalar in a parametric form, R{r), 
from its definition via the metric, given by Eq. (|14p . 
Then, using a specific equation of state pr = Pr{p) or 
Pt = Ptif), one may in principle obtain F{r) from the 
gravitational field equations. Finally, once T = T^^i^i is 
known as a function of r, one may in principle obtain 
f{R) as a function of R from Eq. ([3|). 



J 



A. Traceless stress-energy tensor 

An interesting equation of state is that of the trace- 
less stress-energy tensor, which is usually associated to 
the Casimir effect, with a massless field. Note that the 
Casimir effect is sometimes theoretically invoked to pro- 
vide exotic matter to the system considered at hand. 
Thus, taking into account the traceless stress-energy ten- 
sor, T = —p + pr + 2pt — 0, provides the following differ- 
ential equation 



F" 



b\ b'r + b-2r 
^ r) 2^ 



F' 



2r3 



-F = 0. 



(24) 



In principle, one may deduce F{r) by imposing a specific 
shape function, and inverting Eq. (I14p . i.e., R{r), to find 
r(i?), the specific form of f{R) may be found from the 
trace equation 

For instance, consider the specific shape function given 
by b{r) — r^/r. Thus, Eq. (|24p provides the following 
solution 



F{r) = Cisinh 
+C2 cosh 



V2 arctan 



V2 arctan 



.(25) 



The stress-energy tensor profile threading the worm- 
hole is given by the following relationships 



p{r) 
Pr{r) 



Pt{r) 



Ci sinh 



^/2 arctan 



ro 



\Jr'^ - rl 



2C2A/2(r2 - r2) + 3roCi sinh 



+ (^2Ciy^2(r2-r2) + 3roC2^ cosh 



r 



(C2V2(r2-r2)+roCi)sinh 



+ [c^^2{r'^-rl)^r^C^ cosh 




(26) 



(27) 



(28) 



One may now impose that the above stress-energy ten- considering the values Ci = and C2 = —1. 
sor satisfies the WEC, which is depicted in Fig. [fl by ^p^^jg^ ^j^^p^ i^^c\:^o^ considered above, the 
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FIG. 1: The stress-energy tensor satisfying the WEC, for 
the specific case of the traceless stress-energy tensor equa- 
tion of state, and for the values Ci = and C2 = —1. We 
have considered the dimensionless quantities WECl — rop, 
WEC2 = ro{p + pr) and x = r/ro. 



FIG. 2: The specific form of f{R), for the specific case of the 
traceless stress-energy tensor equation of state, by imposing 
the values Ci = and C2 ~ —1- The range is given by 
< R/Ro < 1. 



Ricci scalar, Eq. provides R — ~2r'^/r'^ and is now 

readily inverted to give r = (— 2rQ/i?)^/'*. It is also con- 
venient to define the Ricci scalar at the throat, and its 



inverse provides ro = (— 2/i?r 



Nl/2 



Substituting these re- 



lationships into the consistency equation ([31), the specific 
form /(i?) is given by 



In principle, as mentioned above one may deduce Fir) 
by imposing a specific shape function, and inverting Eq. 
(fTi)) . i.e., i?(r), to find r{R), the specific form of /(i?) 
may be found from the trace equation ([3]) . In the follow- 
ing analysis we consider several interesting shape func- 
tions usually applied in the literature. 



f{R) = -R{ Ci sinh 



arctan 



1 



+C2 cosh 



arctan 



1 



,(29) 



which is depicted in the Fig. [21 by imposing the values 
Ci = and C2 = -1. 



B. Specific equation of state: pt = ap 

Many of the equations of state considered in the litera- 
ture involving the radial pressure and the energy density, 
such as the linear equation of state Pr = ap, provide very 
complex differential equations, so that it is extremely dif- 
ficult to find exact solutions. This is due to the presence 
of the term F" in pr- Indeed, even considering isotropic 
pressures does not provide an exact solution. Now, things 
are simplified if one considers an equation of state relat- 
ing the tangential pressure and the energy density, so that 
the radial pressure is determined through Eq. (|17[) . For 
instance, consider the equation of state pt = ap, which 
provides the following differential equation: 

P'(^-l)-i^ib-b'r{l + 2a)]^0. (30) 



1. Specific shape function: 6(r) — r^/r 



First, we consider the case of b{r) 
30l) yields the following solution 



rp/r, so that Eq. 



(31) 



The gravitational field equations, p^ - p^ . provide 
the stress-energy tensor threading the wormhole, given 
by the following relationships 



Prir) 



X [2{r^ - rl) + 3ar^ - 4r^a - r^a^] , (32) 



Pt{r) = ap{r) 



Cir'^a 



(33) 



One may now impose that the above stress-energy ten- 
sor satisfies the WEC, which is depicted in Fig. [3l by 
imposing the values Ci = — 1 and a = —1. 

As in the previous case of the traceless stress-energy 
tensor, the Ricci scalar, Eq. ([SI), is given by i? = 
— 2rg/r'^ and its inverse provides r = (— 2rg/i?)^/'*. The 
inverse of the Ricci scalar evaluated at the throat inverse 
is given by rg = (— 2/i?o)^^^. Substituting these relation- 
ships into the consistency relationship Q, provides the 
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1.2 1.4 1.6 1.8 2 2.2 2.4 2.6 2.8 3 



FIG. 3: The stress-energy tensor satisfies the WEC, for the 
specific case of the equation of state pt = ap and considering 
the form function b{r) — r^/r. We have imposed the val- 
ues Ci = — 1 and a — —1, and considered the dimensionless 
quantities WECl = rgp, WEC2 = rl{p + pr) and x = r/ro. 



f{R)/ROo2: 




R/RO' 



FIG. 4: The profile of f{R) is depicted for the specific case 
of the equation of state pt = ap and considering the form 
function b{r) = rg/r. The values Ci = — 1 and a = —1 have 
been imposed, with the range given by < R/Ro < 1. 



specific form of f{R), which is given by 



fiR) 




(a^ + 2a + 2) + (a + 2) 



(34) 



This function is depicted in Fig. [5] as J{R)/Rq as a func- 
tion as R/Rq, for the values Ci = — 1 and a ~ —1. 

2. Specific shape function: b — y/riyr 



Consider now the case of 5 = 
yields the following solution 



/r^, so that Eq. ([30] 



(35) 



The stress-energy tensor profile threading the worm- 
hole is given by the following relationships 



pt{r) = ap{r) 



Cia (1 - y/? 
2r2 



Prir) 



Giro 
' 16r3 



1 - 



(36) 



10J— + ^I^ (14a + 10) + (4a2 - 26a + 5) 



.(37) 



Rather than consider plots of the WEC as before, we 
note that it is possible to impose various specific values 
of Ci and a that do indeed satisfy the WEC. 

Following the recipe prescribed above, the Ricci scalar 
is given by i? = ^Jr^jr^l"^ and is readily inverted to pro- 
vide r = (y^frQ I R)"^ I ^ . The inverse of the Ricci scalar 
at the throat provides ro — I/a/Rq. Substituting these 
relationships into the consistency relationship ([2]), the 
specific form /(i?) is finally given by 



f{R) 



Ci 



i?i -2(RR^Y +R§ 



-L 1 \ 2 ' 

R^ - i?5 



3_9ci -21 + lOa 



8RR^ + ( 11 + 10a ) + ( 2 - 22a + ia^ )r^RS + (-5 + 12a - 4a^) R^ Rq 



(38) 



3. Specific shape function: b{r) — ro + 7^ro(l — ro/r) tion 

17^~2q-1 

Finally, it is also of interest to consider the specific Fi^) — Ci (r — 7 rp) ^ ^^^^ 

shape function given by b{r) — rg +^'^rQ(\ — ro/r), with ^ -,^(i+2a)-i 

< 7 < 1, so that Eq. ([^ provides the following solu- _ j.^)^ (39) 
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It is useful to write the last equation in the form F{r) = 
CiX"r-("+i)y", where X, Y, u, v are defined as 

X r - 7^ro , y = r - ro , 

_72-2a-f _72(1 + 2q!)-1 
"~ 2(72 -f) ' 2(72 -f) • 
I 



As in the previous example, we will not depict the plot of 
the functions, but simply note in passing that one may 
impose specific values for the constants a and Ci in order 
to satisfy the WEC. 

The Ricci scalar, Eq. pi)) , provides R — 27^rQ/r'* and 
is now readily inverted to give r = (27^rQ/i?)^/"'. The 



Thus, the stress energy tensor profile threading the 
wormhole is given by the following expressions: 



(40) 



(41) 

I 

Ricci scalar at the throat is given by Rq = 27'^/rQ, and 
its inverse provides rg = 7 ■\/2/ _Ro ■ Substituting these 
relationships into the consistency relationship ([3]), the 
specific form f{R) is given by 



Prir) 



9l 

2r3 



+r-(2+")^2^2(l0a2+4) 

+r"(i+"V^72u(4Q + 6) 
+2r2-"u (w + 1) 
-r^(i+"V^72u(4a-6) 



" (2a2 + 6a + 4) + r^^^+^Vo (-7a + 2a^-f'^ - 87^ - 7a72 - 2a^ - 3) 
r-"rov (-7^(5 + 4a) - a - 5) + Ar^-'^v (1 + a) 
2r""ro7^w [v - a) + 2r^^°'rQV {-v + 7^ + 1) 
r-"row (4a + 5) (7^ + l) - '^r^^" {u - a) 
+ X'^-^Y'" [2r~'^rl-i^u [u - I) + 2r^-"roU (1 - u) + 2r^-°'u {u - 1) 



J^U—l^V 



u—l\rv — l 



4r "ro7^ww + 4r^ "rouv (7^ + l) - 4r^ "uw 



fiR) 



CiR 



{RoRY 



(a + l) 



7'-(t)M7^ + l) + (f 



rI 



1 -f^-2Q-l 



1 -f-^(l-2a)-l 



27^(a2 + 2a + 2) - ( :|^y (3a + 4)(7^ + l) + ] (2a + 4) 



(42) 



IV. SUMMARY AND DISCUSSION 



In general relativity, the NEC violation is a fundamen- 
tal ingredient of static traversable wormholes. Despite 
this fact, it was shown that for time-dependent worm- 
hole solutions the null energy condition and the weak 
energy condition can be avoided in certain regions and 
for specific intervals of time at the throat [31|. Nev- 
ertheless, in certain alternative theories to general rel- 



ativity, taking into account the modified Einstein field 
equation, one may impose in principle that the stress en- 
ergy tensor threading the wormhole satisfies the NEC. 
However, the latter is necessarily violated by an effective 
total stress energy tensor. This is the case, for instance, 
in braneworld wormhole solutions, where the matter con- 
fined on the brane satisfies the energy conditions, and it 
is the local high-energy bulk effects and nonlocal correc- 
tions from the Weyl curvature in the bulk that induce 
a NEC violating signature on the brane Another 
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particularly interesting example is in the context of the 
D-dimensional Einstein-Gauss-Bonnet theory of gravita- 
tion where it was shown that the weak energy con- 
dition can be satisfied depending on the parameters of 
the theory. 

In this work, we have explored the possibility that 
wormholes be supported by f{R) modified theories of 
gravity. We imposed that the matter threading the 
wormhole satisfies the energy conditions, and it is the 
higher order curvature derivative terms, that may be in- 
terpreted as a gravitational fluid, that support these non- 
standard wormhole geometries, fundamentally different 
from their counterparts in general relativity. In the anal- 
ysis outlined above, we considered a constant redshift 
function, which simplified the calculations considerably, 
yet provides interesting enough exact solutions. One may 
also generalize the results of this paper by considering 
$' ^ 0, although the field equations become forth order 
differential equations, and become quite intractable. The 
strategy adopted to solve the field equations was essen- 
tially to specify 6(r), and considering specific equation 
of state, the function F{r) was deduced from the grav- 
itational field equations, while the curvature scalar in a 
parametric form, R(r), was obtained from its definition 
via the metric. Then, deducing T — T'^^ as a function 
of r, exact solutions of f{R) as a function of R from the 
trace equation were found. 

Furthermore, we note that f{R) modified theories of 
gravity are equivalent to a Brans-Dicke theory with a 
coupling parameter w — 0, and a specific potential re- 
lated to the function f{R) and its derivative. In this 
context, it was shown that static wormhole solutions 
in the vacuum Brans-Dicke theory only exist in a nar- 
row interval of the coupling parameter [2ll |. namely, 
—3/2 < w < —4/3. However, we point out that this re- 
sult is only valid for vacuum solutions and for a specific 
choice of an integration constant of the field equations 
given by C{w) — —l/{w + 2). The latter relationship was 
derived on the basis of a post-Newtonian weak field ap- 
proximation, and it is important to emphasize that there 
is no reason for it to hold in the presence of compact 



objects with strong gravitational fields. 

Another issue that needs to be mentioned, is that the 
above-mentioned interval imposed on w was obtained by 
considering the violation of the WEC (recall that the 
WEC imposes p > and p + Pr > 0) . Now the authors 
in [21] obtained the respective constraints on w by con- 
sidering negative energy densities, i.e., p < 0. This is 
not a necessary condition, as one may consider positive 
energy densities and in alternative impose the condition 
p + Pr < 0, which violates the WEC, and consequently 
the NEC. Note that this is justified as the fundamen- 
tal ingredient in wormhole physics is the violation of the 
NEC, and not the imposition of negative energy densities. 
In principle, this condition combined with an adequate 
choice of C (w) could provide a different viability and less 
restrictive interval (including the value w = 0) from the 
case of —3/2 < w < —4/3 considered in [2l|. 

For the vacuum case considered in the present paper, 
we note that there are no viable solutions, as now we 
have three gravitational field equations and two arbi- 
trary functions, b{r) and F{r), so that the system is over- 
determined. This difficulty can be surpassed by consider- 
ing the general case of $'(r) ^ 0, but now it is impossible 
to find an exact analytical solution, and numerical meth- 
ods are needed to solve the system of equations. However, 
in the presence of matter things are totally different, as 
this adds additional degrees of freedom. Thus, in princi- 
ple one may construct a whole plethora of wormhole solu- 
tions (a specific equation of state was considered in 2'^), 
in addition to adequately choosing C{w) in Brans-Dicke 
theory. Work along these lines in presently underway. 
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